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Abstract 

The Ising approximation of the Heisenberg model in a strong magnetic 
field, with two-, and three-spin exchange interactions are studied on a Husimi 
lattice. This model can be considered as an approximation of the third layer 
of ^He absorbed on the surface of graphite (kagome lattice). Using dynamic 
approach we have found exact recursion relation for the partition function. 
For different values of exchange parameters and temperature the diagrams of 
magnetization are plotted and showed that magnetization properties of the 
model vary from ferromagnetic to antiferromagnetic depending from the value 
of model parameters. For antiferromagnetic case magnetization plateau at 1/3 
of saturation field is obtained. Lyapunov exponent for recursion relation are 
considered and showed absents of bifurcation points in thermodynamic limit. 
The Yang-Lee zeros are analyzed in terms of neutral fixed points and showed 
that Yang-Lee zeros of the model are located on the arcs of the circle with 
the radius R = 1. 

1 Introduction 



Since the spin-1/2 ^He atoms have week attractive potential and light mass the 
zero-point fluctuations are large. Therefore helium can solidify only at very high 
pressures. Unlike usual solids, where the dominant is dipolar magnetic nuclear 
interaction in solid ^He nucleon-nucleon interaction is dominant, therefore nuclear 
spins become ordered at 1 mK (for usual solids 1 microKelvin) [1] . For such systems 
theory of magnetism is based on multiple-spin exchange mechanism. 

It is important to examine solid and fluid ^He films absorbed on the surface 
of graphite, [H IH EJ [6] since it is a typical example of a two-dimensional frustrated 
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quantum-spin system [71 [S]. The first and second layers of tliat system form trian- 
gular lattice, while the third one forms a system of quantum 1/2 spins on kagome 
lattice[9l HOl E]. Both experimental [121 1131 Hi] and theoretical [IH M\ studies in- 
dicate that three site exchange interaction is dominant in these systems. When 
the density of ^He nuclei decreases the magnetic properties of system changes from 
ferromagnetic to antiferromagnetic one. This behavior can be explained in terms 
of multiple-spin exchange (MSE). For high densities the system is fully packed and 
three-site exchange interaction is dominant, therefore the system is ferromagnetic. 
When density decreases two-site exchange interaction becomes dominant and mag- 
netic properties of the system change into antiferromagnetic one. 

The magnetic properties were studied in Bose-Einstein condensates of ultracold 
atoms. The dynamical creation of fractionalized halt-quantum vortices in Bose- 
Einstein condensate of sodium atoms have been demonstrated pTT]. Moreover, there 
were made the description of the Josephson effect in Abelian and non-Abelian Bose- 
Einstien condensate of alkali elements [1 8 [ [T9] . Last of them may be realised in 
superfiuid ^He Josephson weak link[20]. The Kosterlitz-Thouless transition with 
random magnetic field [21] was used to describe marginal metal to insulator transition 
in disordered graphene[22]. 

At low temperatures antiferromagnetic quantum atoms of ^He can exhibit plateaus 
on the magnetization curve. Beginning at some value of the external magnetic field 
(less than saturation field) the magnetization of the system does not change when 
the external magnetic field increase and system adsorbs energy without any change 
of magnetization. When the external magnetic field reaches a certain value the 
magnetization changes its value again. 

This phenomena was theoretically predicted by Hida[23j for a ferromagnetic- 
ferromagnetic-antiferromagnetic Heisenberg chain 3CuCl2-2 dioxane compound, which 
consists of the antiferromagnetic coupled trimers. For the different models the mag- 
netization plateaus are predicted to occur in chains, ladders, within the dynamical 
and transfer matrix approaches (see Ref. [21] -[36]). For two dimensional systems 
appearance of magnetization plateaus has been found experimentally. On the trian- 
gular lattice a magnetization plateau was observed at m/m^ = 1/3 for compounds 
like C6Eu[STlEH], CsCuCls [M] (see also Ref. [10]- [IS]). For kagome lattice also have 
been observed magnetization plateau at 1/3 of saturation field [13 HH]. Noticeable, 
that plateaus appear only for the values of magnetization which are quantized to 
fraction values of the saturation magnetization. This phenomena was theoretically 
explained by Oshikawa, Yamanaka and Affleck[49j in 1997. 

As mentioned above the third layer of the ^He films absorbed on the surface of 
graphite is kagome lattice. Usually, the antiferromagnetic kagome lattice is investi- 
gated using numerical simulation|5U t 1 ^ 15^ . In this paper dynamic system approach 
based on exact recursive relation for the partition function have been used. The key 
point of the approach is the so-called recursive lattice approximation [53] . which is 
a very powerful tool in investigating many theoretical problems in statistical me- 
chanics. We approximate kagome lattice by Husimi lattice which is a Bethe-type 
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recursive lattice. This approach allows us to plot magnetization curves for dif- 
ferent values of temperature and exchange parameters. For the finite size lattice 
at low temperature there are bifurcation points near plateaus, therefore Lyapunov 
exponent [53] |55| |56| |57] of recursion relation are studied to verify disappearance of 
the bifurcation points in thermodynamic limit. 

In the strong external magnetic field Heisenberg model can be approximated by 
Ising one. There are not reasonable conditions in solid and fluid ^He under which 
we could neglect the non-diagonal Heisenberg interaction, nevertheless, if the strong 
magnetic field is directed along the z-axis, we expect that it reduces the transverse 
fluctuations. It is supposed that in this case and cr^-spin components are infinitely 
small and these spin components can be neglected[Tl EZ]- 

The knowledge of the partition function is very important in statistical mechan- 
ics, since the thermodynamic functions of the system can be expressed by means of 
the partition function. When the system undergoes phase transition, some thermo- 
dynamic functions (such as free energy) become nonanalytic at that point, therefore 
phase transitions can be associated with zeroes of partition function. In 1952 for 
the first time Yang and Lee [SHI El] offer a method for studying phase transitions by 
mean of the partition function zeroes. The new concept of the partition function 
zeroes on the complex magnetic field plane was introduced ( Yang-Lee zeros). They 
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studied the partition function of the Ising model as a polynomial in activity (e ''s^ ^ 
where H-is the complex magnetic field) and proved a circle theorem which states 
that for Ising model the zeroes of the partition function lays on the unit circle on 
the complex activity plan. It was shown that in thermodynamic limit the system 
undergoes phase transition only when distribution of zeros on the complex activity 
plan cuts the real axes. 

The Yang-Lee zeros can be studied using the dynamical systems approach jGOj 
EH |62] or transfer matrix method |63| [63] |65| 166] . According to Biskup et a/. [67] 
and Monroe [SB] the Yang-Lee zeroes correspond to the phase coexistence lines on 
the complex magnetic field plane. Phase coexistence lines correspond to the points 
where recursive function absolute derivatives in two fixed points are attractive and 
equal [61]. 

This paper is organized in the following way. Section 2 is devoted to the investi- 
gation of the two-, and three-site exchange interactions Heisenberg model on kagome 
lattice in an external magnetic field. In Section 3, for the strong external magnetic 
field Heisenberg model is approximated by Ising-like one on Husimi lattice and ex- 
act recursion relation for the partition function and magnetization are derived. In 
Section 4 the magnetization curves for different values of temperature and exchange 
parameters have been plotted. The absents of bifurcation points at low temperature 
have been shown by studying Lyapunov exponent. In Section 5 Yang-Lee zeroes are 
studied using dynamic system approach. Finally, Section 6 contains the concluding 
remarks. 
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Figure 1: a) kagome lattice b) Husimi lattice 

2 Two-, and Three-Site Exchange Interaction Heisen- 
berg Model for Fluid and Solid ^He on Kagome 
Lattice 

The Hamiltonian for ^He consists of two parts 

H = Hex + Hz, (1) 

where Hex is spin exchange interaction Hamiltonian and Hz is the Zeeman Hamil- 
tonian which is responsible for magnetism. In the most general form[l] the multiple 
spin exchanges Hamiltonian can be written as 

Hex = -Y.^na{-lYPn, (2) 

where the summation runs over all permutations of particles, Pn is the permutation 
operator of n particles, J„a is the corresponding exchange energy (a distinguishes 
topologically inequivalent cycles), and p is the parity as defined in permutation 
group theory, i.e., it is odd (even) if the decomposition of the permutation into a 
product of pair transpositions involves an odd (even) number of transpositions. For 
third layer of ^He absorbed on the surface of graphite (see Fig. 1 (a)) with two- 
and three-site exchange interactions Hex takes the following form 

Hex = ^2 ^ Pij-^-i E {Piik + P^jk)^ (3) 
Pairs Triangles 

where Pij is a pair transposition operator, Pijj, is an operator that makes cyclic 
rearrangement in the triangle and first sum goes over all bonds while second sum 
goes over all triangles. The expression of pair transposition operator P^ has been 
given by Dirac[2] 

P^J = ^ (1 + '^''^^■) ' (4) 
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one can write 



where cTj is the Pauh matrix, acting on the spin at the i-th site. Using Eq.(|l]) one 
can derive expression for Pij^ 

1 1 

Pijk = PijPik = ^ (1 + ^i^j) (1 + ^i^k) = ^(1 + CTiCTj + CTiCTfc + (cTj CTj ) (cTjCrfe) ) , 
P^jl = PikPij = ^ (1 + ^i^k) (1 + CTjCTj) = ^(1 + O-iO-j + (Ti(Tk + {(T i(T k) {(T i(T j)) . 

(5) 

Using identities 

{CT i(T j) ■ {(T i(T k) = (7 j(Tk + (Ti[(T j X (Tk] , 

{aiau) ■ {(TiCTj) = (Tk(Tj - (Ti [a-j X cr^] , (6) 

Inserting equations (jl]) and ([7j) into ([3]) one can obtain the exchange Hamiltonian 
with two- and three-site exchange interactions 

Hex = yY1(^ + CTiCTj) - y (1 + '^i'^i + ^j^k + (Tk(Ti)- (8) 

The expression for the Zeeman Hamiltonian is 

Hz = -T.l^Ba„ (9) 

i ^ 

where 7 - is the gyromagnetic ratio for ^He nucleus, B - is the magnetic field. 
Final expression for two - and three-site exchange interaction Hamiltonian is 

if = y ^ (1 + (TiCTj) - y XI (1 + ^i^i + ^O^k + CTfeCTi) " ^ ^hB(Ti. (10) 

(ij) {hj,k) i 

3 Ising Approximation of Heisenberg Model on 
Husimi Lattice. Recursion Relation for The 
Partition Function 

Several approximations can be applied in the Heisenberg Hamiltonian (fTOj) . At first 
the Pauli matrices in ( JTOl) can be replaced by classical three-dimensional vectors s 
of unit length (so called 0(3) classical Heisenberg model). Moreover, in the strong 
external magnetic field aligned along the z-axis the contribution from x and y com- 
ponents of classical spin variables will be insignificant and the main contribution 
will be from the z component, which can effectively take values = ±1, so instead 



5 




Figure 2: The procedure of derivation of the recursion relation for the partition 
function on Husimi lattice 

of Heisenberg model we have the Ising one. Consequently instead of Hamiltonian 
fflOl) we get, 

if = y ^ (1 + SiSj) - y (1 + ^i^j + ^j^k + SkSi) - h ^ Si, (11) 

where h = ^KB^. 

The kagome lattice can be approximated by Husimi lattice which is a recursive 
one (see Fig. 1). The recursive lattice gives an opportunity to obtain the exact 
recursion relation for the partition function and apply dynamic system theory. The 
recursive Husimi lattice is formed in the following way: to each site of the first 
(central) triangle 7 — 1 triangles are attached. The construction of Husimi lattice 
continues recursively for each triangle. If 7 = 2 then Husimi lattice is an approxi- 
mation of the kagome one. (see Fig. 1 (b)). 

Using the hierarchical structure of Husimi lattice one can derive recursion relation 
for the partition function. The partition function of the system with Hamiltonian 
(UnD is 

where ks is Boltzmann constant, T is temperature and the sum goes over all 
{si, S2, ...Sn} configurations of the system. In the rest of the paper all constants 
will be taken in Boltzmann constant's scaling {^2/^3 = J2-, Js/ks = J3, h/ks = h). 
Using Eq. f lTT]) one can write 

H {Si, S2, ...Sn) "^2 /. , \ I "^3 , , , \ , ^ 
-TPf 2^ (1 + SiSj) + — 2^ (1 + SiSj + SjSk + SkSi) + 7f2^S,. 



kbT 2T 2T,^^, ' ' - rj. 



(13) 
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To obtain recursion relation for the partition function one can separate Husimi lattice 
into three identical parts (branches) and at first realize summation over all spin 
configurations on each branch, then to sum over spins of the central triangle. The 
result of the summation for each branch will only depend on value of corresponding 
spin variable (see Fig. 2). By denoting 

A(si, S2, S3) = -y (3 + S1S2 + S2S3 + S3S1) + y (1 + S1S2 + S2S3 + S3S1) , (14) 
we can rewrite the expression for the partition function in the following form: 

7 \^ I ^ l^Sl j ^l^Sl +^2 +S3 M (0). , (0)x . (0)^ 



(15) 

where sum goes over central triangle, Qnisf^) denotes the contribution of a branch 
at the Si'^^-th site of the central triangle and n is number of generations. 

At the same way one can express Qnisf') in terms of gn-i{sf'') cutting it along 
any site of the first generation. Therefore 

g^is^)- E exJ ^\ ' ^ \ gn-As^')9n-As^\ 

(16) 

After the summation we get the following expressions for gn{si) 

/ \ -3J2+2J3+2fe 9 , , ^ _£2 , , , , -J2-2h 

(17) 

where we denote gn{si) by 5'n(=t) depending on Si sign. By introducing Xn — 
9n{+)/gn{—) variable the recursion relation for the partition function can be ob- 
tained 



1 + 2ca: + ac^x'^ 
^^'^ ^ a + 2cx + c^x^ ^ ^'^^ 

where a = e * c = e^. The thermodynamic functions of the system, such as 

magnetization, can be expressed in terms of Xn- The magnetization of the sf'^ site 
is expressed 

^ = (^S°^> = e-^is,,s2...s.)/ksT ■ (19) 

{Si} 



Let us at first realize the summation over central vertex Si and then over all other 
spin variables. After inserting equations f lT7|) into f lT9|) the formula for magnetization 
will take the following form 



For arbitrary value of magnetic filed h, with given temperature and exchange param- 
eters one can draw the dependence of magnetization from external magnetic field 
by implementing the simple iteration from the recursion relation for f{x), begin- 
ning with some initial value of xq, The thermodynamic limit correspond to infinite 
number of iterations (n — )■ oo). 



4 The Magnetization and Lyapunov Exponent 



In this section the magnetic properties of the model in an external magnetic field 
have been studied. As known, J2 and J3 are not obtainable in the experimental 
measurements, because n-spin exchange also makes a contribution to {n — l)-spin 
exchanges, but there is some effective exchange parameter, J = J2 — 2 J3 which can 
be directly obtained from the experiments. The magnetic properties of the model 
depend on value of effective exchange parameter and vary from ferromagnetic into 
antiferromagnetic one, depending on whether two- or three-spin exchange interaction 
is dominant. From experimental measurements and theoretical calculations we know 
that at low-density region three-site exchange interaction on the regular triangular 
lattice is dominant. The corresponding estimated value ^1] of exchange parameter 
is J = J2 — 2J3 = —3.07mK. According to general principles ground state of the 
system in this case is ferromagnetic. Magnetization function of the model with 
exchange parameters J2 = 2 rriK and J3 = 2.5 mK is presented in Fig. 3. For 
high temperatures the magnetization curve has monotone form of Langevin type 
(Fig. 3(a)). The value of saturation field as large as higher temperature (Fig. 3(b)). 
For some value temperature (critical temperature) the magnetization curve ceases 
to be smooth, the jump of magnetization takes place at an arbitrary low value of 
the applied magnetic field (Fig. 3(c)). The ground state of the model ordered 
ferromagnetic. The magnetization in zero field is equal to its maximal value, which 
corresponds to the ferromagnetic phase with all spins pointed in the same direction. 
Under the effect of an arbitrary weak magnetic field the twofold degeneracy of this 
phase is removed by orienting all spins along the field. The corresponding diagram 
is presented in Fig. 3(d). 

The model is antiferromagnetic if two-site exchange interaction is dominant. In 
Fig. 4 are presented magnetization curves that correspond to the values of exchange 
parameters J2 = 3 mK and J3 = 2 mK. As expected the magnetization curve is of 
Langevin type for the high temperatures (Fig. 4(a)) and the saturation field decrease 
with the temperature (Fig. 4(b)). With further decreasing of the temperature 



m = 




(20) 



8 




1.0: 
0.5 ; 
B 0.0: 

-0.5 ; 

-1.0- 



c) 



r 



-0.5 



0.0 

h 



0.5 



-1.0 



1.0 



1.0 
0.5 

B 0.0 

-0.5 
-1.0 

d) 



-0.5 



-1.0 -0.5 



0.0 

h 



0.0 

h 



0.5 



0.5 



1.0 



1.0 



Figure 3: The magnetization curves for J2 = 2mK, J3 = 2.5'mK, n = 50000 and for 
(a) T = 1.5mA", (b) T = 0.7mK, (c) T = O.GmK, (d) T = 0.3mK 



plateau at 1/3 of saturation field appears on the magnetization curve (Fig. 4(c)). 
This phenomena can be explained in the following way. For each triangle two spins 
directed along magnetic field and the other one oriented opposite to the field ("up- 
up-down" phase) . 

At lower temperatures the magnetization plateau turns into bifurcation points 
and period doubling (Fig. 4(d)). These bifurcation points will disappear if number 
of iterations tends to infinity. As it is known Lyapunov exponent |54[ [55| [56| |57] 
is exactly zero simultaneously with bifurcation points, indicated in some external 
parameters and temperature. The Lyapunov exponent is the index of exponential 
divergence of two near points after N iterations: 



ee 



N\{xo) 



/(^)(xo + e)-/^"^(xo) 



(21) 



N 



where f^^\x) = f{f{. ..fix))). In the limit at e ^ and ^ cx) ([2lD gives the 
exact formula for X{xq) 



Afxn) = lim lim — In 



fW(xo + e)-f(^\xo) 



lim — In 

N-foo N 



df^'^Kxo) 



dXn 



■ (22) 
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Figure 4: The magnetization curves for J2 = 3mK, J3 = 2mK, N = 50000 and for 
(a) T = ImK, (b) T = O.SmK, (c) T = O.lmK, (d) T = O.QbmK 



Using formula for derivative of composite function 

= fifixo))f'{xo) = /'(xi)/'(xo), Xi = /(xo), 



XO 



one can transform Eq. (!22|) in the following form 



A(xo) = lim — In 



_d_ 

dXr 



f 



lim — In 



N-1 

n /'(^^ 

i=0 



(23) 



(24) 



If Lyapunov exponent is negative, then the final state of the system after iterations 
tends to stable point of stable cycle composed of more than one points. The positive 
Lyapunov exponent result in chaotic final state of the system and zero Lyapunov 
exponent correspond to bifurcation points. Using Eq. one can plot dependence 
of Lyapunov exponent from external magnetic field. In Fig. 5(a) Lyapunov exponent 
curve for J2 = 3 mK, J3 = 2 mK and T = O.OdmK has been plotted. Near plateau 
Lyapunov exponent is very close to zero but it does not cut the axes. In Fig. 5(b) 
the magnetization curve for the same values of exchange parameters, temperature 
and number of iterations n = 1 500 000 000 has been plotted. As can be seen from 
figure the bifurcation points disappeared. 
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Figure 5: The plots of Lyapunov exponent (a) and magnetization curve (b) for the 
values of exchange parameters J2 = 3mK, J3 = 2mK, temperature T = 0.05mK 
and iterations n=l 500 000 000. 

5 Yang-Lee Zeroes On Husimi lattice 

The thermodynamic properties of the system may be investigated by studying the 
dynamics of the corresponding recursive function (|T8|) . As known from dynamical 
system theory when the number of iterations tends to infinity {n — 00) the resulting 
point can have different behavior: 1) limiting point tends to one point, 2) limiting 
point tends to the set of points 3) limiting point has chaotic behavior. In the first 
case resulting point can tend to one of the stable fixed points. The point x* called 
fixed point of recursion relation Xn = if x* = f{x*). The fixed point can 

be stable (the iterations of any point xq near x* tend to the x*), or not stable (the 
iterations of any point xq near x* tend to move away from x*). The stability of fixed 
point depends on value of the derivative of f{x) at that fixed point x*. If |A| < 1 
where A = f'{x*) then the fixed point called attracting, if |A| > 1 then the fixed 
point called repelling and if |A| = 1 the fixed point called neutral (indifferent). The 
values of external parameters (temperature, magnetic field, etc.) at which system 
has only one attracting fixed point correspond to a stable paramagnetic state. If 
the system has two attracting fixed points then the stable state corresponds to the 
fixed point with maximum value of |A|, moreover the values of external parameters 
at which |Ai| = IA2I correspond to two possible ferromagnetic states with opposite 
magnetizations. 

The fixed points of the recursion relation ( !T8|) determined by following equation 



X 



1 -|- 2cx + ac^x"^ 
a + 2cx + c^x^ 



(25) 



In general this equation has three complex solutions. If there are only one attracting 
fixed point then the system state is paramagnetic and phase transition does not 
present. The values of external parameters at which system has two attracting fixed 
points correspond to metastable region|6^. The boundary of metastable region can 
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Figure 6: The border of metastable region (dashed curve) and phase coexistence 
hues (sohd curve) for J2 = 2mK, J3 = 2.5mK and (a) T = 1.5mK, (b) T = 0.9mK, 
(c) T = O.QmK, (d) T = O.bmK, (e) T = O.SmiT, (f) T = O.lmK. 

be found from the condition that one of the fixed points become neutraL There are 
not phase transition on the boundary of metastable region. The phase coexistence 
hues can be found from the condition that absolute values of derivatives at the fixed 
points become equal[6^. When phase coexistence line cuts the real axis then the 
system undergoes first order phase transition. The locus of the phase coexistence line 
inside metastable region correspond to the locus of Yang-Lee zeroes. Consequently, 
the boundary of metastable region can be determined from the following system of 
equations: 

f{x) = e*^ • ^ ' 

Eliminating x from this equations for given temperature and exchange parameters 
one can find the equation for c parameter. 

In Fig. 6 by dashed line boundary of metastable region for J = —3mK (J2 = 
2mK, J3 = 2.5mK) and different values of temperature are plotted. By solid line 
phase coexistence lines inside the metastable region are plotted . If T > Tc phase 
coexistence line does not cut the real axes and is an arc of a circle with radius R = 1 
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Figure 7: The border of metastable region for J2 = 3mK, J3 = 2mK and (a) 
T = ImK, (b) T = 0.5mK, T = O.lmK. 

(Fig. 6 (a), (b)). For T < Tc phase coexistence hne cuts the real axes and the first 
order phase transitions occurs, therefore system has ferromagnetic behavior (Fig. 6 
(c), (d), (e), (f)). Our calculations show that for this values of exchange parameters 
Te = 1.2426 ■ ( J3 - J2) = 0, Q213mK. 

If the two-site exchange interaction is dominant the system has antiferromagnetic 
behavior. For this case in the metastable region there are two attracting fixed points 
but they never become equal (no phase coexistence line in the metastable region), 
therefore Yang-Lee zeroes correspond to the boundary of metastable region. In 
Figure 7 boundary of metastable region calculated using equation (!26l) for J2 = 
3mK, J3 = 2mK and various temperatures are plotted. There is a relation between 
appearance of plateaus and locus of the Yang-Lee zeroes. The locus of the Yang- 
Lee zeroes (Fig. 7 (a),(b)) for the values of temperature for which plateau does 
not appear on magnetization curve (Fig. 4(a),(b)) have been bounded. The locus of 
Yang-Lee zeroes is unbounded (Fig. 7(c)) if the plateau appear on the magnetization 
curve (Fig. 4(c)). 

6 Conclusions 

In the present paper dynamic system theory has been used to study solid and fluid 
^He films absorbed on the surface of graphite. The third layer of ^He films, which 
is kagome lattice was approximated by (7 = 2) Husimi one. In the strong external 
magnetic field the Ising model have been considered instead of Heisenberg one with 
two- and three-site exchange interactions. This approach allows us to obtain mag- 
netization curves for different temperatures. The magnetization plateau at 1/3 of 
saturation field has been obtained for some values of exchange parameters and tem- 
peratures. By studying Lyapunov exponent for antiferromagnetic case the absence 
of the bifurcation points on magnetization plateau at low temperature have been 
shown. 
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The Yang-Lee zeroes have been studied in terms of neutral fixed points of the 
recursion relation. It was shown that Yang-Lee zeroes on the complex e^'^/"^ plane 
located inside the metastable region (existence of two attracting fixed points) and 
correspond to phase coexistence lines (the lines where absolute derivatives in fixed 
points are equal). The locus of the Yang- Lee zeroes for different values of exchange 
parameters and temperature have been plotted. It was shown that for the ferromag- 
netic case Yang-Lee zeroes are located on the arc of the unit circle. For antiferro- 
magnetic case if there are magnetization plateaus the locus of Yang-Lee zeroes are 
infinite, otherwise they are finite. 
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